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ABSTRACT 


In the present work, a study of the optimum design 
of R.C. Intze tank is attempted. The entire superstructure 
and suhstructur e of the water tower is divided into three 
parts viz. the tank, the supporting frame and the foundation 
and its optimum design is obtained interactively by solving 
three nonlinear programming problems sequentially » 

Two problems viz. 600 KL oap a3ity~''tank with 15 m 
staging and 1000 KD capacity tank with 22m staging are solved 
and the number of columns and the number of bracings at 
the optimum for each problem are obtained by parametric 
study. 
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CENTER 1 


&ENERAI( INTRODUCTION 


1,1 Introduction: 

Outstanding progress Das been achieved in the field 
of devising better design solutions, in the last two decades 
at the advent of the electronic digital computer. 1 conscious 
engineer has to investigate several alternative designs and 
has to choose the best among them and it is not a 30 b of 
great difficulty at this age of computer. Many possibilities 
obviously exist, with two clearly defined extremes. One 
extreme is to utilize the computer capability to the fullest 
and ’automate' this search; the other extreme is to utilize 
human intution in an 'interactive manner' to guide the 
computer in its calculations. Structural optimization is a 
broad interdisciplinary field which requires skillful 
combining of mathematical knowledge with engineering and 
it IS both intellectually attractive and technologically 
rewarding. Hence any attempt to exploit the field of 
optimization for the development of better designs is 
worth investigating. 

Extensive literature is available in the field of 
structural optimization as a result of simultaneous aeuelopment 
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towards rigorous structural analysis, and efficient 
algorithms for optimizations. Package programs for the 
automated optimal design of various structures like 
Transmission line towers, bridges, multistorey frames have 
been developed and are in use. Structural optimization study 
for chimneys, silos, bridge decks. Transmission towers and 
so many others structures has been made and encouraging results 
have been observed. In this w&rk an attempt is made to 
investigate the optimum design of reinforced concrete Intze 
type water tank, 

1.2 Present Literature: 

Though there are numerous papers available for the 
approximate and rigorous analysis for various types of tanks » 
the literature is practically nil towards the minimum cost 
design of water tank structures. But there are some general 
indications in selecting the type of the tanks from the various 
available types ( 6 ) like square or nectangular tanks, 
cylindrical tanks (with slab bottom or with spherical dome 
bottom), Intze tank etc, O.P, Jain (pf) has indicated that 
the Intze type tank is the economical one in the capacity 
range between 225 KL and 1800 KL. O.P. Jain and Jai Krishna (8) 
have indicated that the ratio of the diameters of the cylinder 
and bottom ring beam should be around 0,7 and the ratio of the 
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Cylinder diameter and the total height of the container 
should he around 1*0 for the economic design of the container, 
ipproximate analysis of the staging frame for wind load (0) 
and for earthquake (1^) are available* Charles Anthony 
111 by (ife) has given the state-of- the art regarding the analysis 
of various types of tanks (rectangular, cylindrical with 
or without monolithic connections) based on elastic theory 
and plastic theory. M.M. Basel e and IT.L. Khusalini (2. ) have 
developed equation to predict the cos"^ range for cylindrical 
t anks , 

1.3 Scope of the Present ■'!?ork; 

In this work computer programs are developed to get the 
minimum cost design of Intze type water tank* The study is 
motivated because for a given capa,city and height of staging 
under the specified environment the number of columns and 
number of bracings which leads to the minimum cost design is 
not known and so also the geometry of the container which 
gives the minimum cost of the container. Furthermore, the 
question arises that after selecting the type of the foundation 
what IS its configuration to get the minimum cost. The present 
study addresses itself to answering the questions raised above. 

Chapter 2 describes the analysis of the container 
both for membrane and effects of continuity, the analysis 



of the staging frame both for earthquake and. mnd forces 
and the analysis of the raft slab foundation. 

Chapter 3 gives the detailed formulation of the optimal 
design problems for the container, the supporting frame and 
the foundation. It describes the objective function (cost 
of the structure) and the reasons for neglecting some variables 
from the design vectors. The behaviour constraints and the 
geometric constraints on the optimal design and a method to 
coalesce the bounded constraints are explained therein. 

Chapter 4 gives the details of the minimization method 
used, A brief description of the sequential unconstrained 
minimization technique (STMT), the procedure for unconstrained 
minimization and the one dimensional minimization method used 
IS given in this chapter. The algorithm used to predict the 
initial interval of unccrtainity is also explained therein. 

The limitations and potentialities of the developed programs 
are indicated. The necessary modifications required for 
the use of the program in some akin situations are hinted. 

Chapter 5 gives the details of the numerical output 
of this work and the conclusions based on the study of the 
results of the developed programs. 



CHIPTER 2 


MiLYSIS MD DESIGN OE WATER TOWER 

2.1 Introduction: 

The complete water tower structure can be considered 
as a combination of the three substructures viz,, (i) container 
or tank, (ii) staging or supporting frame and (iii) foundation. 
The method of analysis of these substructures and the forces 
and moments required for the design of the various elements 
of the substructures are discussed in brief in this chapter, 

2.2 Analysis and Design of the Container; 

The Intze tank consists of four shells i.e, two 
spherical? one conical and one cylindrical shell which are 
connected by three ring beams as shown in Fig^ 1, The loads 
and forces on each shell are symmetrical about the axis of 
revolution and the pure membrane state of stress will exist 
so long as each shell is simply supported at its edges i.e, 
it is able to undergo resulting edge displacements without 
restraint, while the supports supply the necessary reaction 
to balance the meridional forces. But in the Intze tank, at 
the junction of shells, the edge displacements are actually 
restrained. Thus the complete analysis requires both membrane 




BB-RING BEAMS 


IG1 GEOMETRIC AND STEEL VARIABLES OF THE CONTAINER 
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and continuity analysis. This is well docuin«bed in (8) 
and a hrief description follows for compfsteness of 
presentation, 

2,2.1 Membrane analysis: ^ ^ 

(1) '1’op spherical shell: The top spherical dome can be 
designed for the meridional thrust at the edges and for 
the maximum hoop force. 

Meridional thrust T^ = wR / ( 1+cos 0 ) 

2 

Maximum hoop force = wR(cos 9 + cos 9--1)/(1+cos 9) at 9=0 
where w = L.L + D.l, 6 is the semi central angle of the dome an<3- R 
Is the radius of the dome. 

(2) Ring Beam 1: The meridional thrust, T^ , exerts a vertical 
load T^ sin 9 on the wall and also it imposes an outward 
radial force T^ cos 9 which causes hoop tension in the 

ring beam 1 . 

Hoop tension = T-| cos 0 D/2 
where, D is the diameter of the tank, 

(3) Cylindrical tank wall: The wall is assumed free to 
deform at both edges and is thus under a pure hoop tension. 

The maximum hoop tension can be found at the bottom and the 
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wall Gan be designed to resist this hoop tension. As the 
pressure wanes linearly with depth, hoop tension will also 
rery linearly. 

Hoop tension at any depth = p^ D/2 , where p^ is the 
pressure at depth d. 

(4) Ring beam 2; If the weight transmitted through the tank 
wall at top of conical dome is w per unit length and 0 is the 
inclination of the dome with the vertical, then the meridional 
thrust caused is w.sec 0 per unit length of circumference. 

Outward radial force in ring beam 2= sin 0 = w tan 0 

Therefore, Hoop tension in ring beam 2 = w.tan 0,I)/2 
The ring beam should be designed for the hoop tension 
calculated as above. 

( 5) Conical dome : The conical dome supports a uniform 
vertical load from walls at its top edge. At top of the 

dome a hoop tension will be created which exerts a radial inward 
force at each slanting strip (dome is assumed as consisting 
of individual slanting strips monolithically joined) and will 
oppose its rotation outward. The magnitude of the radial 
force created at top edge is so much that on combining with 
the vertical load, the resultant lies along the meridian 
of the conical dome. Thus the vertical load at top edge of 
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tl e conical dome is supported by it with the creation of 
m jridional thrust and a hoop tension. The water pressure on 
conical dome and its own weight acting at any point give rise 
to hoop tensions at each plane, whose inward reaction, together 
with the water pressure and weight of dome, cause a resultant 
force which is meridional. 

Meridional thrust at any height d = (% D'cos 0) 

where is the total vertical load acting on the dome above 

that height and 1’ is the diameter at that level. These 
expressions are used to calculate the necessary forces to 
design the sections of the conical dome, 

(6) Bottom dome: like the top sperical dome, the bottom dome 
also develops only compressive stresses both meridionally and 
along hoops. The maximum meridional thrust at the edges and 
maximum hoop force (at 0 = 0) can be calculated as described 
earlier. 

(7) Ring Beam 3 : Ihis beam receives an inward inclined thrust 
from the conical dome and inclined outward thrust due to 

the reaction from the bottom sperical dome. The vertical 
components of these two thrusts add up and the ring beam is 
designed for this load. Their horizontal components oppose 
each other and depending upon their relative magnitudes, ring 
beam is either in hoop tension or in hoop compression and can 
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"be designed accordingly, 

(q) 

2.2,2 Continuity analysis:^ 

Tlie forces due to continuity are obtained by applying 
tbe principle of consistent deformations. The vertical displace- 
ments are always consistent at each joint as each is shell 
is free to defonn in this direction and consistency has only 
to be satisfied for horizontal and angular displacements 
between the shells meeting at a 30 int. Based on the stiffness 
of each shell at edges for horizontal and angular moments, 
equations of consistency are framed and the forces due to 
continuity can be determined after solving this system of 
simultaneous equations in rotation (6) and displacement (x). 


Membrane deformations and stiffnesses at the edges of 
various elements can be e3q)ressed as follo\/^®: 


(l) Top dome: Slope at left edge 2wR sin 0/Et( clockwise) 

2 

wR sin 0 

Horizontal deflection x, = — (VCl+cos e)-Cos 6) (Inwards) 

^ Et 

RBt 1 

Moment stiffness per unit length — (k + r” ) 

^ 4 V ^ 


2 

Corresponding radial force per unit length = Et/(2„Mc sin 0) 

2 

Thrust stiffness per unit length = Et/(^Rlc sin 9) 

2 

Corresponding moment per unit length M 2 =Et/C 2 >k: sin 0) 
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4^2 1 

where , ^ = 3 (“t ) ; k = 1 - cot 9 

^ 2 A 

(2) Ring beam 1: Radial thrust to cause unit outward 

E hd 

deflection — 


R 


Moment per unit circumference to cause unit rotation 

Ebd5 


M3 = 


12 R 


(3) Oylindrical wall: Moment stiffness = 2^z 


Corresponding thrust Ih = 2/o. z 


■4 = r 

3 

Thrust stiffness = 4 ^ z 
Corresponding moment M3 = 2 ytl 


pR 


Displacement of tank wall at "bottom x, = 

w Et 

PR.^ 


Slope of the wall = 

^ Efch 


(4) Ring "beam 2; Radial tlirust to cause unit outward 

Ebd 


deflection ^ 

R 


Ebd- 


Moment per unit circumference to cause unit rotation M^r = 5- 

12R^ 


where, b and d are breadth and depth of the beam and R is the 
radius of the cylindrical shell. 


(5) Conical dome : Outward deflection at the top edge — 


D-ELfc 


2 Et 


Outward deflection at the bottom edge = 


2 IJt 



12 


tan 9 

Slope at the edges = (2T2-T^) 

Stiffnesses at the top edge: 

Et.k, 1 

Moment stiffness M, 


(k^ .k^-k2k^) tan^e 


Corresponding thrust = 


Et k2 


( k^k^-k 2 k^) sin 0 tan 0 


Thrust stiffness 


% 


Et k. 


1 


1 sin e (k^k^-k 2 k^) 


Corresponding moment Mg = 


Et. k. 


sin 9 tan 9(k^k^-k2k^) 
Stiffnesses at the bottom edge: 

Moment stiffness 


Et k^ 1 


Corresponding thrust, Hg = 


(k^ k^ -k2k^)tan 0 
Et k^ 


(k^k^-k^kp sin 0 tan 0 


Thrust stiffness H. 


Et. k» 

__L 


10 p 

I'sin e(k'k’ -k^k') 
14 "^3 


Corresponding moment M^q = 


Et,k> 

0 


(k’k'-k'k') sin 0 tan 0 
^14 2 3 


For the given geometry of the conical shell k^,k 2 ,k^,k^ and 
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’ ^2 ’ ^4 obtained from fhe Table 1 shown, 

after evaluating 8 ,8’ as below. 1 is the slant length 

of the cone (complete) and 1* is the truncated cone slant 
T ,, 4 12 cot^e 

length. ^ = - — — — where 6 is the semi apex angle of 

t-^ 

the cone. 


e = 2Afl and 8’ = 2^^!’ 

(6) Ring beam 3: Ingular rotation 0 


12 R d 

(m-H — - ) 


bd^B 


where, m is the clooicwise edge moment per unit circumference 

length applied to the beam at its top and H is the outward 

radial force, b and d are the breadth and depth of the ring beam3. 

d 

Inward moment of top edge due to an angular rotation = 0 — 


Outward moment due to hoop tension = 


H R 
bd E 


R 


Equating these two, 0 = 


Moment stiffness , = 


5R 


bd^ E 
bd^E 


3R 


m 


per radian 


Corresponding outv/ard thrust = 


bd^E 


2R‘^ 


per radian 


bdE 


Thrust stiffness H^2~ 


per unit movanent 


R" 


Corresponding edge moment M ^2 “ " 
(7) Bottom dome: Slope at edges 


bd^E 


2R‘ 


per unit movement 


2 w R sin e w R sin 0 

5 = - -£ 


Et 


Et 
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TABLE 1 : OOElCil, SHELL STIPENESS GOISTIHTS 


^ or e* 


^4 


hb 

^4 

10 

183.0 

25.10 

6.79 

169.0 

23.55 

7.30 

11 

220.5 

27.25 

7.21 

236.5 

29.35 

7.90 

12 

288.0 

32.60 

7.91 

312.0 

36.00 

8.70 

15 

366,5 

38.40 

8.60 

396.0 

42.20 

9.40 

U 

463.0 

46,20 

9.36 

490.0 

48.85 

10.09 

15 

568,5 

52.50 

10,10 

596.0 

56.30 

10.82 

16 

696.0 

60.30 

10.82 

728.0 

63.20 

11.42 

17 

834.0 

68.40 

11.52 

880.0 

72.30 

12.13 

18 

990.0 

76,65 

12.24 

1038.0 

81.20 

12.97 

19 

1172.0 

86.00 

12.98 

1222.0 

89.70 

13.61 

20 

1370.0 

95.90 

13.69 

1430.0 

100.00 

14.30 

21 

1574.0 

105.60 

14.39 

1658,0 

111 .0 

15.10 

22 

1830.0 

116.3 

15.10 

1900,0 

121.1 

15.70 

23 

2100.0 

127.8 

15.88 

2190.0 

134.0 

16.50 

24 

2400.0 

139.9 

16.45 

2470.0 

145.0 

17.16 

25 

2695 

175.0 

17.24 

2800.0 

156.9 

17,80 
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2 

wS sin 0 1 

Horizontal moment of edges x = - — >( - — _ cos e) 

1+cos 0 

_2 3 

w R sin 0 w R sin 0 

_ ^ 2cos20+cos0-3) 

2 Et 6 Et(1+cos 0) 

RJH: 1 

Moment stiffness M. , — -(k + r — ) 

4 ?^ ^ 

2 

Corresponding thrust = Et/(2;\k sin0) 

Thrust stiffness = lt/(2>^k sinO) 

2 4-2 

Corresponding moment = Et/(2Ak sin0) where, )v=3(R/t) and 
1 

k = 1 cote . 

2> 

Reactions due to continuity at ring beam 1 junction: 

Let the net rotation of the joint be 5^ clockwise and 
net displacement be inward. Then the changes in the slope 
and displacement from the membrane states can be expressed 
as follows: 


EL ement 

Slope 

Displacement 

Top dome 

5l - Sd 

X- Xd 

Ring beam 


Xi 

Tank wall 

«r«w 



where, x^, x^ displacements of dome and wall respectively and 
6dJ 6^ slopes of dome and wall respectively. 
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Iluations of consistency of deformations (joint at Ring iDeam 1) : 

I^(x^-x^) + + H^(S^-6^+H^x^+( thrust from dome)=0 

These two equations can he solved for 6-| and and the reactions 
due to continuity can be evaluated, to R^ are the reactions 
derived from these evolution. 


Element 

Moment 

Thrust 

Hoop tension 

Dome 

^1 



Ring beam 1 

®2 

^*5 

% 

Wall 


% 

'>9 


R^ = M^(6^-6^) + E2=M3d^ ; R5=M4(6^-6^ + H^x^ 

~ B[2^^1'”^d^ ^ 5 R6'~^4^ 

x.B JL, x.Eig x^E A 3 

where , A^t A 3 are the cross sectional area of the respective 
member for unit length . 

Reactions due to continuity at the ring beam 2 joint: 

let 62 and X 2 be the rotation and displacement of 
this joint, then changes in the slope and deformation of the 



17 


elements and their reactions can he expressed as follov/s; 


EL ement 

Slope 

Displacement 

Moment 

Thrust 

Hoop Tension 

Wall 


"2-^w 

^10 

R ^5 

^16 

Ring 

heam2 

^2 

X2 

^11 

^4 

R17 

Conical 

dome 

62-6^ 

^2-^0 

R-,2 

®15 

^18 


6^, - slope and displacement of the wall 

displacement of the conical dome 


+ H^(x:2-x^4-Mg52+M^(x2-x^)+Moment due to halcony slah=0 

M5(62-6^+H5(D^-x^+H6Xg+l^(62-6c)^-I^(^-^P = 0 

These two equations can he solved for 62 
reactions can he evaluated as shown helow. 


R^3 Mj( 62 - 6 ^+H^(x 2 -x^ ; R^^-HgX2 ; 62 “ 6 Q)+Hg(x 2 -x^) 



(P/2) 



. 

(D/ 2 ) 



^ 2 ^ 

(D/2) 


where A^,4^,lg are the cross sectional areas of the respective 
elements for unit length. 
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Re ctions due to continuity at the ring iDeam 3 30int: 

let the rotation and displacement of this ^oint be 5^ 
and then the changes in slope and displacement of the elements 
(from the membrane state) and their reactions can be expressed 
as follows. 


El ement 

Slope 

Displacement 

Moment 

Thirust 

Hoop teusion 

Gonical Dome 

^3"^c 

X3-XQ 


R22 

R25 

Ring beam 3 

63 

2:3 

^0 

^3 

^6 

Bottom dome 

63-6^ 

X3-X13 


^24 

^7 


6-^ and Xi^ represents slope and displacement of the bottom dome 
respectively. 

Hio(6^“6q) + H-j >|(x^— x^)+M^2S^+H^2^3 

+ ( Thrust from both the domes) = 0 

These equations can be solved for x^ and 6^ and the reactions 
to R^rj can be calculated, 

R-jg “* ’ ^0~^1 1 ^ 3'^^1 1 ^3 

R 2-1 “ M^g(6g— 6-^) + g( Xg— X|^) ; R22“^-|0^ 1 ^^3"”^c^ 
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R 25 » ^26 ^7 calculated using the corresponding 

X*? SA. 

unit length areas in as "before, where Do is the diameter 

(V2) 

Of the ring beam 3. 


The continuity alters the hoop forces near joints and 
imposes meridional moments and they die out after some distance 
and hence they are of local nature maintaining the meridional 
stresses in the interior portions unaltered. The design will 
be improved by providing sufficient reinforcements and 
thichnening the shell edges, based on the above calculated 
reactions. 


Analysis of the ring beam3: 

The load of the tank is transferred to the columns 
through this circular beam. This beam is suboected to a 
uniformly distributed load all along its length. The 
following expressions are used to calculate the design 
moments of the circular beam. 

2 

Support moment = P wR (2a) 

2 

Mid-span moment=P ' wR ( 2a) 

2 

Maximum tv/isting moment =P’ 'wR (2a) 


where 2a is the angle subtended at the centre by two column 
points (ie.) 2 a = (m/n^^) radians, vhere nQ-number of columns* 
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p = (1-a c Ota) /2a ; =(-1 .0+asinc(+cos^a/sina)/2a 

p’'= (/_o^-acos0-acota a )/2a , where a =sin^a-cosa(a^^sin^a) V2 

o S 

and 0 = sin \ a ) 
s 

P,p' , p’ ’ are the curved heam coefficients which depends on 
the number of columns (i.e,) a at the centre. 

2,3 Supporting frame inalysis: 

Idle supporting frame consists of R. 0, columns, 
placed symmetrically on the circumference of a circle and braced 
together by horizontal members as shown in fig, 2, The columns 
are subjected to the vertical load and the lateral forces 
due to wind or earthquakes, acting on the tank and the frame. 

The problem of analysis can be done rigorously by considering 
it as a space frame but in the normal design procedure and in the 
situations where repeated analysis is required one should go 
in for the ^proximate method of analysis, 

C«) 

2.3*1 Ifi-nd load analysis: 

Since the tank base, foundation and the braces are 
very stiff compared to the columns it is considered that 
the tower deflects, maintaining the column axes almost vertical 
at their top, bottom and at their junctions with braces 
developing the points of inflexion at mid heights of each 




FIG-2 GEOMETRIC AND STEEL VARIABLES OF THE FRAME 



FIG 2 b SHEAR AND MOMENT IN THE 
BRACING BEAM 
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panel, Tlie internal reactions at these inflexion points can 
he calculated as described helow. The deflected pattern of 
the tower is shown in Pig, 2(a), 


Considering the total cross sectional area of all the 
columns as the area of circular ring, moment of inertia Qf 2 

the assumed cantilever beam of circular ring section =■ 


8 


8 


Therefore, vertical force on the farthest leeward side 


column = 


4m 

n^a I) 


a = 


4m 

n^D 


(Because 


m p 


n^aD^ 2 
8 


= stress ) 


where, m — moment in the assumed cantilever beam, nj^-number 
of columns, a- area of one column, D-diameter, Vertical force 
in columns lying on bending axis is zero. 


Shear stress in the column(at 9 from the bending axis) 


Q 


q. = 


2t sec0 n^ a I) /S 


2 2Q cos 0 

2(l)/2) t cos9 = 




Shear force maximum '( at 0 = 0) S-q a 


2Q 


, where 


n. 


IS the lateral force acting at that level. 


The maximum moment will occur at the top and bottom 

h 

end of each panel and will be equal to S ~ where h is the 
clear panel height. Thus these forces and moments (i.e.) 
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vortical load, bending moment and shear force can be used for i 

the design of columns, 

[ 

Analysis of Bracing beams: , 

The shear forces acting on either side and the 
moments set up at the joint are shown in the Big. 2(^1.) , Considering ' 

f 

the static equilibrium of moments at the joint, 

m-] nip M 

= is the equation 

sin^e+it/nQ) sin(0-'n:/nQ) sin(27t/nQ) 

C Q^h^+Qoho) 2 

relating m^ , nip and M. But = — cos 9 whereQ^ , Qp 

and h-^ , hp denote shear forces and panel heights above and 
below the bracing, respectively. 


Design forces and moments for the bracing beam: 

Maximum moment = — — sin( 9^^^Ti/ng) where, 

sin( 2%/nQ) 

^^ax " M ® “ ®max* 


9 = can be obtained from the condition for maximum 

HI SIX. 

moment m^(i,e,) 9 which satisfies 2 tan9 tan(9+'n;/nQ)=1, 

4 * 2 2 

Shear force maximum = (cos 9 sin(0+'Tt/nQ)-cos (9-2Tc/nQ) 

MB n„sin(2'n:/n^) 


sin (9-3n:/nQ)) 

where, 0 — and AB is the length of the bracing beam. 
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■j'he twisting moment can be assumed as 5 percent of the 
iraximum bending moment. 

2.3.2 Earthquake load analysisr 

The supporting tower should be designed for the forces 
and momenis due to earthquake loading also. The design forces 
can be obtained by any one of the two methods specified in 
IS=rl893 (criteria for earthquake resistant design of structures) 
i.e. 1, Seismic coefficient method. Response spectrum method. 

In the later one the earthquake factor is calculated and the 
earthquake force is obtained which is coming as shear in the 
columns, 

Sa 

Earthquake factor a, “ 

n o g 

p = coefficient depending upon the soil- foundation system 
I = coefficient depending upon the importance of the structure 
seismic zone factor for average acceleration spectra as 
given in Table 2 of 13:1893. 

Sa 

average acceleration coefficient depending on the damping 
and appropriate natural period, 

I =1.5, P = 1.0, are the coefficients considered in this 

work and 5 percent damping is considered as per 5.2,2 IS: 189 5 
specification for elevated tower supported tanks, 

Natural period of vibration T= sec. 
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roM 

whore , -A== f dx : 

El 


E - Youngs modulus of concrete 
I - Moment of inertia of the column 


Es rthquake force “ * (Total load +1/3 Self load) 


Shear force per col-umn = (Earthquake force) /Number of columns 


0 

Moment in the column = sheao* force* —k— where h.-clear panel heigh 
(Meq) ^ 

Reaction force in the column due to the moment at that 

M 

level can be obtained by P = — r R and then the direct load on 

the column is added to it to get the total load coming on 
the column* The column should be designed for these two forces 
and P^^. The moment in the bracing beam can be calculated 

^ M. ^4. 

from and the beam can be designed accordingly. 


Moment in the bracing beam ^br 



cos(ix/2-'n:/nQ) 


the number of columns. 


where n_ denotes 
c 


2,4 Foundation Analysis and Design: 

Depending on the soil characteristics the selection of 
the type foundation can be made. There are various types of 
foundation which can be provided for water tank stmeture, 

A combined raft foundation with complete raft slab and the 
foundation with ring raft slab(i.e,) slab combined with the 
bottom circular beam, are the two types among them. Pile foundation 
is also sometimes used. The analysis of ring type of raft 
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as shown in Fig. 3 which is selected in the present work is 
■briefly explained below: 

The uniformly distributed load coming on the circular 
beam is calculated. Then the design moments (i.e.) support 
moment, midspan moment and twisting moment can be calculated 
using the circular beam coefficients P , P* and p'* as shown 
earlier. The raft slab can be considered as a cantilever slab 
with the soil reaction as the tiniformly distributed loading. 

The maximum negative moment can be calculated and the slab can 
be designed for these forces. The reinforcements can be 
curtailed at the half distance keeping the alternative bars to 
the complete length, and the distribution steel can be provided. 
The actual pressure of the soil should be less than the 
allowable bearing capacity of the soil. This analysis is an 
approximate one and is q^uite suitable for embedment in the 
optimal design routine because if is simple and reasonably 
accurate and will not consume more time during the function 
and gradient evaluation, 

2,5 Summary: 

The membrane analysis and continuity analysis of 
the tank, approximate method of analysis of the staging frame 
for wind forces and earthquake forces, and the analysis of 
raft slab foundation have been explained in this chapter, 

T hese are embedded in the optimization programs explained 
in the succeeding chapter. 
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OPOTiiL DiSIGN FORHULATIOF 

3*1 Introduction: 

The structural optimization problem can be fomulated 
either using mathematical programming approach or the optimality 
criteria approach* The present problem of optimizing the 
design for minimum cost of a water tower structure is 
formulated as a mathematical programming problem* The design 
vectors, constraints and the formulation of the objective 
functions of each problem are esplained in this chapter. 

The one-shot fomnulation of the optimum design of 
the entire water tank stmcture becomes a complex problem 
due to the following facts. 

1, Since the optimum solution seeking is a repetitive analysis- 
design procedure the analysis of the composite structure viz, 

the analysis of the container with membrane and continuity 
effects, the frame analysis of the supporting tower for wind 
and earthquake loadings and the foundation analysis is time- 
consuming procedure. 

2, Since the entire water tank structure consists 
of numerous variables, the problem becomes a large 
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non-linear programming problem which requires an efficient 
solution technique, sufficient computer memory and more 
computer time for seeking the solution. 

3. Since the number of columns and the number of bracings 
influence the entire analysis-design procedure, as well as 
the size of the design vector and the problem formulation, 

it is difficult to include them in an automated optimum 
design procedure, 

4. Since the type of the foundation differs for different 
soil conditions it is difficult to formulate the problem 
including the foundation unless the analysis-design of the 
various foundations are developed with the proper strategy to 
select the foundation. 

Therefore, in the present work optimum design of the 
water tank is obtained by substructuring the problem into 
three parts viz, the tank, the staging and the foundation 
and solving them sequentially. 

5,2 Optimum Design Formulation of the Container: 

The geometry of the container of the intze tank is 
shown in Fig.1, t^ ,t 2 ,t^,t^,t^ are the thickness of the top 
dome, wall thickness at top edge, wall thickness at bottom 
edge, conical dome thickness and the bottom dome thickness 
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1 aspect ively.D^ and are the diameters as shown in Pig, 1. 
" ^ , d^, d 2 , h^, d^ are "breadths and depths of the 

lUng "beam 1, 2, 5 respectively and h^,h 2 ,h^,h^ are the 
heights and rise of the domes as shown in the Pig, 1 . 


The various steel reinforcements are also shown . 

Thus there are fifteen steel variables and seventeen geometric 
variables in the problem. To reduce the complexity of the 
problem the design vector is considered to have only fifteen 
variables. The reasons for neglecting the remaining variables 
from the design vector are e3cplained below: 


(1) The top dome carries its own weight and the 
nominal live load and in general the design comes out with 

the minimum requirements of steel and thickness as per IS code, 

(2) The geometry (breadth and depth) of the top 
ring beam can be assessed judiciously which depends on the 
hoop tension on the ring beam from the top dome. The reinforce- 
ments Can be provided for this hoop tension based on the 
allowable stress in steel. 


(3) In order to reduce one variable due to the 

variation In the oylinarioal wall, the bottom edge thiotaiess 

* -ui o -r.-pn< 5 f'ribine: minimum thickness 
13 oonsidered as a design variable prescribing mi 

based on code reauirement at the top edge. 

(4) Curved beam oooffioients shows that the maximum 

positive banding moment at the mid span is approximately 
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V. qual to half of the support moment. Since the cross sections 
of the circular beam at the support section and at the mid 
span section are same, half of the reinforcements provided 
for the midspan moment will be bent and placed at top 
of the beam from which it can be seen that the beam will 
have equal area of steel both at top and bottom at the 
support section. This relation between them facilitates 
to remove one steel variable from the design vector. 

(5) The rise of the top dome and the bottom dome 
can be fixed as 1/8th of the C 7 linder diameter and 1/6th 
of the bottom dome diameter to avoid h^ and h^ from the 
design vector. 

(6) The reinforcements in the bottom dome generally 
comes as per minimum requirement and hence they are neglected 
from the design vector, 

(7) The steel variables like the distribution steel, 
stirrups can be assessed properly based on the code requirement, 
and hence they are not considered as design variables. 

Thus the above variables which seems to have negligible 
effect on the cost of the structure are neglected. The design 
variables considered in the problem are listed below (Ref,E'ig,1) . 
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3.2,1 Design variables: 


Diameter of the tank (D^) 

height of the cylindrical shell (h 2 ) 

thiclaiess of the cylindrical shell at the bottom edge(t^) 

breadth of the ring beam 2 (b^) 

depth of the ring beam 2 (dg) 

thickness of the conical shell (t^) 

height of the conical shell (h^) 

diameter of the conical dome (D 2 ) 

breadth of the ring beam 3 (b^) 

depth of the Ring beam 5 (d^ 

thickness of the bottom dome (t^) 

area of hoop steel at the cylindrical shell bottom edge 

(Ast^) 

area of hoop steel at the Ring beam 2 (Ast^) 
area of hoop steel at the conical dome (ist^) 
area of steel at the support section of the ring 
beam 3 ( Ast^-|) , 


The area of steel in each case represents the fraction 
of respective cross-section of the elements. 


3.2,2 Oboective function: 

The objective function is the cost of the tank which is the 
function of the design variables, The total cost of the tank 



53 


<■ an be expressed as the sum of the cost of concrete, steel 
and shuttering. 

Total cost fCx!) = Cost of concrete 0^+ Cost of steel 

+ Cost of shuttering 

*5 

where, C^, C^, denotes cost of concrete per m , cost 

2 

of steel per tonne and cost of shuttering per m respectively. 

Tq, Wg» denotes volume of concrete, weight of steel and 
total shuttering area. Since the cost of shuttering varies 
for conical dome and spherical dome (from the cost for 
cylindrical wall shuttering) proper weight factors have been 
considered while calculating the over all shuttering area 
in terms of straight shuttering area. 

2 

(i.e.) Cost of shuttering for conical dome=aQCgj^ per m 

_ 2 

Cost of shuttering for spherical Pei* m 

= 1.25 and = I .5 are the weight factors considered 

in the present work 

2 Di“I>2 2 

h 2 l'^'tTt(D-|+b^ ) b2d2'^''rt(l'j'^‘3^2'^ ^ 2 ^ ^ "^ 4 . 

1 2 

+ % 2.0. 833. D2 * — - h2.t5+ixD2-'b3«^5 t^/4+7tb^b^d^ 

= the sum of volumes of various elements(i,e.) top dome. 
Ring beam 1, cylindrical wall, Ring beam 2, conical dome, 
Ring beam 5 and conical dome. 
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t.2 

mD 

‘^s ( Ast^4-lst2)4'n:D^b^d^( Ast^-HAst^)+7T:D-| h.^( Ast^+Astg) 

D-1-i-Dp P Du“*D^ p A lo 

+7tD^b2d2(,Ast^+ istg)+ -rtC — ~) t^( Astg+lst^g) 

2 ^ 'I *5 

+71^2 g-0.833t^+bjd^uD2(2lst^^+ ~ Ast^2^ 

p / p 13'i"i“Dp 

A^=27:(D^+t5)h2+a(l(TtD/4i-27xO. 833^2/6.0) +(-~~^)7t “c* 

= the sum of the areas of the shuttering 

3,2,3 Constraints: 

In seeking the optimum design which has the minimum 
value of the above formulated objective function, it has 
to satisfy certain requirements. These requirements turn out 
to be equations and/or inequalities which are functions of the 
design variables. These constraints are of two types. One 
type is geometric or side constraint which can be expressed 
explicitly in terms of design variables. The other type is 
based on the structural response called behaviour constraints and 
which Cannot be normally expressed explicitly. The evaluation 
of the behaviour constraints requires the analysis of the 
structure at the current state. 

The shuttering for conical dome is expensive and if 
the inclination angle is more than 45° one has to provide 
double shuttering, both on top and bottom of the dome. To 
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,,void this, the angle can he provided as less than or equal 

i:o 45°. 

hi 

GJ( 1 ) = = tan 0 < 1 

(D^-D2)/2 

The ratio of the depth of breadth of the ring beam 3 is maintained 
between 1 and 3. This range of dimensions facilitates the 
direct access for the torsional coefficients during analysis 
and it also restricts the depth of the beam to increase to a 
large value. 



This will come as two separate inequality constraints. 


GJ(2) =1.0 - 


h. 


The minimum thickness of the 
should not be less than 10 cm 


0 , GJ(3) = -+ < 0 . 

conical shell and bottom thickness 


GJ(4) = 1.0 - — i < 0; 

10 ■" 


C5J(5) = 1.0 - 


tc 

-SL. < 0 
10 “ 


The minimum thickness of the cylindrical shell at the bottom 
should be more than or equal to 15 can. in upper bound can also 
be prescribed for this variable to avoid very large thickness of 
the wall. 

GJ(6) = 15.0-t5 < 0 ; GJ(7)=t5-100.0 < 0 
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The thickness of the cylindrical wall at the bottom edge 
siould be less than the breadth of the ring beam 2, which 
is an obvious requirement to avoid constructional difficulty. 
The depth of the ring beam should be more than the breadth 
of the beam, 

GJ(8) = 1.0 ~ b2/t5 < 0 , GJ(9) = 1.0 -cL2/b2 < 0 

The geometry of the tank (container) should be such that it can 
store the given volume of v/ater. This requirement can be 
expressed as a strict equality constraint in the optimization 
problem. But to avoid the difficulty which will arise in hand- 
ling the equality constraint in the search procedure, the 
above requirement is expressed as two inqualities as follows. 


The capacity should be more than(or equal to) the 
required volume but it should not be more than 5 percent 


of the required volume. 

GJ(10) = 1.0 


Volume( required) 
Capac ity(provided) 


< 0 


Volume required . 

GJ(11) — - ''.05 < 0 

Capacity provided 


The reinforcements provided in various elements should 
not be less than 0.3 percent and also it should not exceed 

section. This limitation creates 


2,3 percent of the cross 
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t e following set of constraints. 

GJ(12) = 0.003 - St^ < 0 , &J(13) =St^-0.023 < 0 

GJ(14) = 0.003 - St2 < 0 , G-J(15) =St2-0.023 < 0 

GJ(16) = 0.003- Stj < 0 , GJ(17) =St5-0.023 < 0 

GJ(18) = 0.003 -St^ < 0 , GJ(19) =St^-0.023 < 0 

GJ(20) = 0,003 -St^ < 0 » GJ(21) =St^-0,023 < 0 

GJ(22) = 0.003 -Stg < 0 , GJ(23) =Stg-0.023 < 0 

where, St^ is reinforcement in the ring heam 1 , 

St 2 IS reinforcement in the ring team 2 

St^ IS reinforcement in the bottom edge of the cylindrical 

wall 

St^ IS reinforcement in the conical dome 
St^ is steel area in ring beam 3 

Stg is steel area of stirrups in the ring beam 3 per unit 

length, and they are fractions of the design variables, 

GJ(24-) =0,003-x^ 2 » GJ(25) = 0.023 - x^2 

GJ(26) =0.003-x^^ , GJ(27)= 0.023 

GJ(28) =0.003-x^^ , GJ(29)=0.023 - X-,^ 

GJ(30) =0.00 3-x^ 5 , GJ(31)=0.023 - x^^ 

where, x^2»^13»^14»^15 design variables pertaining to steel 
area-' as dataxlGd Barlxsr* 
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stresses in Gorlov.,^+ 

iicrete and steel in fhc 

leabers should be less tb +u 'Various 

eue.es ;:: IC 

can not be ee^ressed explnolty 1 constraints 

Pin ^ function of design variahioo 

However thev cnn -k^ -. ^ variables , 

ey can be evaluated by analysis procedure. 

W(32) = cq/oj . ^ 

SJ(33) =r,/o;. 1.0_< 0 

0 


<W( 34 ) = (^/a^ _ ^ 

ej( 35 ) = o^/o; _ 1,0 ^ 
ej(36) = o^/cr . 1,0 _< 


»J(37) 


^9 

- 1.0 




< 


0 

0 

0 


(JJ( 38 ) = 5 q _ _ Q 

(j_ 

s 


< 0 


GJ( 39 ) = Si 




-“•.0 < 0 


OTp 

W(40) = -11 _i,o 
G-J(41) = -13 , 

GJ( 42 ) = 

0 - 


< 0 


- 1.0 < 0 


- 1.0 


GJ(43) = -. 

( 

ej(44) = -SL _i 0 
0 - 1.0 


°k. 

Sb 


0 
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./here t ^ - allowable direct compressive stress; 50 ksc for 
1200 and 40 ksc for M150, 

0 "^ - allowable tensile stress in concrete; 17 ksc for M200, 

01 - allowable stress in steel: 14 OO ksc, 

allowable compressive stress in bending; 70 ksc for 
M200 and 50 ksc for Ml 50, 

0 "^ - stress in concrete due to meridional thrust in the 
top dome, 

0 2 — stress in concrete due "to hoop force in the top dome, 

or - stress in concrete due to hoop force in the ring beami , 
3 

01 — stress in concrete in the conical dome, 

6 

01 - stress in concrete due to meridional thrust in the 
8 

bottom dome, 

Og - stress in concrete due to bending moment at the 
top dome edge 

O^Q -stress in steel due to bending moment at the 
top dome edge, 

CT^^- stress in concrete due to hoop force at the 
bottom edge of the cylindrical wall, 

(^2~ stress in concrete due to hoop force in the 
ring beam 2 , 

(T^- stress in concrete due to hoop force in the 
1 3 

conical dome, 

stress in concrete( tensile) due to bending moment 
at the edge of bottom dome. 
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stress in concrete due to bending moment at tbe 
bottom edge of conical dome, 

0^6"* stress in steel due to bending moment at the 
bottom edge of conical dome. 

It Can be seen that there are large number of 
constraints which arise due to the lower and upper bounds on the 
design rariable, For example steel area should lie between 
0.003 and 0 , 023 . These kinds of constraints can be modified as 
equivalent constraints which eventually reduces the total 
number of constraints considerably. 


Equivalent constraint for bounded variable constraint: 


Suppose there are two constraints which relate the 
specified values and maximum and minimum of the design 
variable (or any function of design variables), then these two 
constrains can be converted into an equivalent constraint as 
f oil ows , 


as < 4 and d < d^^ 

respectively. 

d^ij^-d 

<0 IS the nomalised form of the first 
max“Siin“" 


d - d^ ^ 

constraint (i.e) \ < 0 and — < 0 is the normalised 




•in 


form of the second constraint. 



4-1 


'"inax"~*\QirL 


\iin^ ^in” ^ax 

wm-l i 

d -d 
max min 


- = «D -1 < 0 

d-d ^ ' 

max min 


In order bo satisfy doth the inequalities, i«e. 

d-d^^ 

< 0 lias to he satisfied . If > 0 is 


n'' n 


^max ^in 


considered as D^, then the equivalent constraint becomes 

V> i 


The above mentioned method of modifying the constraints 
reduces the tot?il number of constrains from 44 to 31 in the 
present problem. 


3.3 Optimum Design Formulation for the Supporting Frame: 

As mentioned earlier the supporting frame consists of 
columns and bracing beams which supports the container and 
transfers the load to the foundation. For a given capacity 
and height of staging one cannot predict easily the number of 
columns and number of panels which gives the minimum cost 
design of the frame. It is usual to provide equal panel height 
though there is no guarantee that the panels of equal height 
will give tho minimum cost design of the frame. But in order 
to formulate the problem as a nonlinear programming problem 
one should predict the number of columns and number of panels 
before hand, otherwise it will have to be formulated as 
an integer non-linear programming problem with the panel 
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heights and number of columns and bracings ais integer variables, 

Sd ace the solution seeking of the nonlinear integer progrannriing 
p oblem IS not easy, the approach of parametric study has been 
adopted in the present work. In other words minimum cost design 
of the staging is studied for various sets of number of 
preassigned columns and braces. Based on the result of this 
study the number of columns and braces corresponding to the 
optimum design are given. 

Suppose ]!„ and IL are the number of columns and number 

c p 

of panels as the preassigned values for a given investigation. 

Then the formulation of the optimum design of staging proceeds 
as follows, 

3,3.1 Design variables: 

The dimensions of the bracing beams and columns and 
the reinforcements provided in these members, are the main 
variables which influences the cost of the supporting frame. 
Diameter of the column (d.) and area of steel (Ast) in the column, 
breadth (b) , depth (d) and the steel areas in the bracing 
beams are bhe variables in the problem. They are shown in 
Big. 2. One can not consider all of them as design variables 
because the three variables in each column, 5 variables in each 
bracing and N variables of panel heights makes the design vector 
a large one. Even for a small value of 1^~3 which is a practical 
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one, the size of design vector turns out to be 22, The design 
\ actor IS reduced considerably by neglecting some variables after 
onsidering the following aspects, 

1. The lateral ties in the columns and the stirmps 

in the beams are not considered as design variables because the 
effect of "fehese variables over the cost is comparatively less 
than the effect of the other steel variables in the column 
and in the bracing beams. 

2. It IS not practicable to go in for the variation 

in the cross-sectional dimensions of the column along its bei^t. 
Therefore, the diameter of the severe column is considered as 
the design variable and this dimension is kept same throughout the 
height of the column, Ttrrther, while the area of steel in 
different panels are considered as design variables, the area 
of steel within a panel is kept constant, 

3. Since the number of panels is preassigned and the 
provision of equal panel height has been accepted in the 
present work , the variables regarding the heights of the panels 
are getting fixed automatically. 

4-. The bending moments in the various level of 
bracings are very nearly equal and they differ only by a 
slight quantity. Thus the dimensions , breadth and depth of 
the lowest level bracing beam can be considered as the design 
variables and the same dimensions can be provided for the other 
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3racing ‘beams, Howerer, to account for tlie variation in the 
bending moment at each bracing level, the area of steel in 
each bracing beam is considered as a design variable. 

5. Depending on the wind direction the nature of the 
bending moment in the bracing beam changes. Hence it is 
conventional to provide equal steel area at the top and bottom 
of the bracing beam. This reduces one more steel variable in 
each bracing beam. 

Thus the design vector consists of the following design 
variables. 

- Diameter of the lower most panel coliomn (d^) 

X 2 - breadth of the lower level bracing beam (b) 

- depth of the lower level bracing beam (d) 

X^, X^. . * steel areas in the columns of each panel , 

Ast^, Ast 2 »Ast^. .. .Istj^p where Ifpy= Np+3 

%p3+1’%p3f2’ %p3+(Np-1) “ steel area in each bracing beam 

of various levels, ( i. e.) ,..Astj^p^2>^s‘^2Hp-1 

The variables corresponding to steel area are expressed as fraction 
of respective cross sections and (Hp— 1) denotes the total number 
of braces. 

It IS to be noted that while the number of columns 
does not affect the size of the design vector, the cost of the 
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-I rame depends on the number of columns provided in the 
t aging, 

5.3.2 Objective function; 

The objective function which represents the cost of 
the staging can be expressed in terms of the design variables 
defined above. 

Total cost of the frame T_ = Vol^ .0^ + W„.C„ 

^ o C So 

7^ o 

where, volume of concrete Vo1q= d^ + (Ilp-I).b.d 

Tl 

(D sin .Nq 

% 

Weight of steel in the columns Wg^ =,7 g( JL d^ -l-( Ast.+Ast^+Ast, 

p 

+ • + Ttdg Ej-.Ng A^s) 

0 TC 

Weight of steel in the bracings Wg-j^= J g(bd(D.sin— "<3 .q) 

^c 

(2Ast(j^p^^^+2Ast(j^p^2)»*-*2Ast(2Up_i) + Stirrup steel)) 

and Wg = WgQ + Vgi^ ; I>= diameter of the ring beam 5; Il|-- height 

of staging; A^-area of lateral ties; s-spacing of the lateral 
ties; Cg-cost of concrete per m ; and C^-cost of steel per tonne. 

3.3.3 Constraints: 

The number of constraints of this problem depends on 
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"the nuniTDer of panel y. Various geometric or side constraints 
and tTae behaviour constraints, on the problem, expressed as 
inequalities are explained below. 


The minimum value of the diameter can be assessed 
depending on the breadth of the ring beam 5 and on the panel 
height in order to keep the column as a short column, 

GJ(1) = 1.0 - < 0 

The breadth of the bracing should not be greater than 
thd diameter of the column and from the practical consideration, 
it should also be slightly lesser than the diameter which may 
facilitate to extend the bars from the beams to the column 
in order to anchor it in the column reinforcements conveniently, 

b < 0,8 d^ (i.e.) &J(2) = 1.0 - 0.8 < 0 


The depth to breadth ratio should be withm 1 and 4 
<i d 

> 1 and 4 are the two inequalities which can be 

D ■” b 

expressed as one equivalent inequality as shown earlier. 


G-J(3) = < 0 where I>^=( 




-1)/(4 -0 


In each panel the columns are subjected to axial force 
and some bending moment. The safe design should satisfy the 
following inequality in each panel column for both wind and 
earthquake loading. 
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(i.e.) 




< 


1 


GJ(4), GrJ(5), ...GJ(lTp) are the constraints of the above form 
each one to the panel colimn and they can not be e 35 )ressed 
explicitLyby the design variables. 


The areas of steel which are expressed as percentage 
of cross section should not be less than the minimnin value 
of 0.008 and it should not be greater than 0,08 (maximum) 

GJ(]Tp+1) = where D^=(Z^-0.008)/(0. 08-0. 008) 

and GJ(Ifp+2), GJ(NpH-3) . . . . .GJ(2lTp-1) can also be expressed as 
above with the corresponding steel variables. Similar 
limitation on the steel variable of the bracing beam forms 
Np-1 number of constraints which are expressed as an enquivalent 
inequivality as shown below, 

G>J(Np2+'l)=\(I^-1.0) < 0 where p^KX^^j^p-0. 00 3)7(0.023-0 .00 3) 

and GJ(Np2+2),GJ(Np2+5),f ....GJ(Np. 5 ) can also be expressed as 
above using the corresponding steel variables IIp 2 = 2 ]Tp- 1 , Np^= 
Np2+-Np-1 and 0.003 and 0.023 are the minimum and maximum steel 
areas, 

lor example if the number of panels Fp=4 , following 
fourteen constraints will be set up. First three for the 
limitations of the diameter of the column, breg-dth of the 



48 


teacmg and the ratio of depth to breadth respectively. The 
next four constraints (Np behaviour constraints) are for 
the stress limitations. The next four constraints (Np) for 
steel area limitations in the columns and the last three 
constrains (Np-1) for the steel area limitations in the three 
(Np— 1) bracings. The total number of constraints are fourteen, 
Seventeen when Np=5, and twenty when Np=6 irrespective of N_. 

5.4 Optimal Design formulation for foundation: 

The type of the foundation is to be selected based 
on various considerations like the soil characteristics and 
the feasibility of construction etc. 1 ring beam-raft 
foundation is considered in this study. The size of this 
problem is comparitively smaller than the container and the 
staging design problem as is seen from the following sections. 

3.4.'! Design variables: 

The foundation consists of a ring beam with raft slab 
and the geometric and steel area variables are shown in fig. 3. 
jSll possible variables are considered as design variables in 
this problem. The design vector consists of nine variables 
which are given below. 

X^- breadth of the ring beam (b^) 

X 2 ‘- depth of the ring beam (d^) 
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X^- thickness of the raft slab at the junction of the 
ring beam (t^) 

X^- thickness of the raft slab at the end (t2) 

X^- projection of the raft slab (12) 

area of steel at the top of the ring beam , Ast^, 
for negative B.M. at the support section 
Irj- area of steel at the bottom of the ring beam, 
for positive B.M« at the middle section) 

2 ^- area of steel at the support section of the 
slab- l 3 t^ 

Xg- area of stirrup reinforcement for unit length of 
the beam, Ast^ 

Xg, Xrjf Iqj are the steel variables which are 
expressed as fraction of cross section of the respective elements. 
Ist^ which represents the distribution steel in the slab is 
not considered as design variable because it can be accounted 
as 0,15 percent of the cross section. 

3 . 4,2 Objective function: 

The cost of the foundation can be expressed in terms 
of the design variables as given below. 

Total cost of the foundation Vol^ C^+Wg 0 ^ where, 

yoi -volume of concrete =tcD(x-| X2+( 2 x^-x^) (x^+x^)/2.0) 

G 



50 


Wg-weight of steel ='n:D( x^X2(xg+0.7x,j)+x^0, 6 Xg4-x^(x^+x^) (3^+0.0015)) 

y - density of steel 
' s 

D - diameter of the ring beam 

3 

C , C - cost of concrete per m , cost of steel per tonne, 
c s 

5.4*3 Constraints: 

This foundation design problem consists of twenty 
inequality constraints out of which 13 ax*e behaviour constraints 
and the rest are side constraints as detailed below • 

The breadth of the beam should be greater than or equal 
to the diameter of the column (dia^). 

b > dia^ ; (i.e.) GJ(1) = 1.0 - Z^/dia^ < 0 

The ratio of breadth to depth of the ring beam should be 

d 

between 1 and 4. These are -^<4- an<i ^ which can be 

b D 

represented as an equivalent constraint, 

1< «i. < 4 ; &J(2) 0 where, Dn=(^ -1.0) /3.0 

The thickness of the raft slab at the end should not be less 
than the minimum thickness of 15 cm. 

t2 > 15 ; GJ(3) = I.O-x/15.0 < 0 

The prooection of the raft slab should not be two large in 
which case it approaches to the complete raft slab type of 
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foundation and the analysis of the raft slab as a cantilever 
may not be valid, 

2d2 < D-dj ; GJ(4) = 2x^/(]>d^)-^1 .0 < 0 whore, 

d^ is the central diameter which checks the closing of the 
two raft slab projections. 

The thickness of the slab t^ should be less than the 
depth of the ring beam 

t^< d^ ; GJ(5) = 1.0 - < 0 

The thickness of the raft slab at the end should be smaller 
or equal to that at the support end 

t^ > t^ ; GJ(6) =^5/x2 - 1.0 < 0 


The actual bearing pressure should be within the safe 
allowable bearing pressure of the soil. This can be expressed 
in terms of the actual provided area and required area to be 
provided to keep the soil with in the allowable limit. 


Actual area Area required ; G<J(7) 1* 



The depth of the ring beam required to resist the bending 
moment as a sin^y reinforced section with the nominal cover 
of say five cms should be greater than (or equal to) the 
arailatle depth and stoilarly with the thiotoiess of the raft 




52 


GJ(8) = 1.0 


depth 


GJ(9) = 1.0 - 


X3 

thictoiess 


The areas of steel e^^ressed as fraction of cross section should 
he within the limits 0,003 and 0.023, 

0.003 < Ast^(Xg)< 0.023 ; GJ(10)= ~ ° 

0.003 < istg (X^) < 0.023; GJ(11)= i ° 

0,003 < ist^ {Iq) < 0.023; GJ(12)= < 0 

0,003 < Ast3 (Ig) < 0.023; &J(13)= *0) < 0 

where, are the steel wariahles as shown in the Jig. 3, 

and = (x^-0,003)/0.020 i=1,2,3,4 

3=6, 7, 8, 9 


The stresses in concrete and steel should he within 
the allowable stress values. The following seven stress limit 
constraints are the behaviour constraints, 

^ S ^ f where, 0^, 0^, 

cr, are allowable stresses in steel, concrete in bending and 
sn 

in tension. 

stress in steel due to the bending moment in the ring beam 
(support section) 

cr^- stress in concrete due to the bending moment in the ring 
beam (support section) 
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stress in steel due to the EM in the ring he^(niidspan section) 

(J^- stress in concrete due to the in the ring beamCmidspan sectic 

CT^— stress in steel due to the BSC at the raft slab support section 

CT^ stress in concrete duQ-to the B*M* at the raft slab support seel 

0^- shear stress in concrete at the support section of bhe slab. 


The above inequalities relating the actual and the allowable 
atresses can be put in the normalized form as follows. 


GJ(14) = 
GJ(16)= 




'5 

■q" 


GJ(18) = 



s 


1*0 < 0 

1.0 < 

1.0 < 0 


0 


GJ(15) = 

; GJ(17) = 




1.0 < 

1.0 < 


; GJ(19) = 

GJ(20) = 


or 


^h 


i-1.0 < 


- 1.0 < 


0 

0 

0 


3.5 Smnmary: 

The formulations of the optimization problems have been 
explained in this chapter. The design variables, the objective 
function and the constraints of each problem is discussed . 

Each one of them turns out to be a nonlinear programming problem. 
These are embedded in the development of the computer program 
which are explained in the subsequent chapter. 
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OPTIMUM SEEKING METHOD iND PROGRiM DEVELOPMENT 

4.1 Introduction: 

Once the problem of optimum design has been formulated, 
the solution of the problem can be obtained using efficient 
algorithms. The optimum design problems as explained in the 
last chapter are nonlinear programming problems. The method 
used for the solution of nonlinear programming problems and 
the details of the developed program are explained in this chapter, 

4.2 Solution Technique: 

There are many algorithms available for the solution of 
nonlinear mathematical programming problem. Penalty function 
method is one among the category of indirect methods. The 
inberior penalty function method is used in this work and the 
well known sequential unconstrained minimization technique (SDMT) 
is briefly described below. 

4-2.1 SUMT algorithm: 

SUMT considers an unconstrained function 0(X,r) where, 

/ V 1 

r IS the penalty parameter, as 0(X, r) = f(x:)-rS_ . 

D-T 

f(x;) Is the objective function and constfaints. 

As r tends to zero ,the minimum of the objective function f(2!) 
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will be the minimum of the unconstrained function 0(X,r), 

The unconstrained minimization is carried out by using 
Davidon-I'letcher-Powell (DFP) variable metric algorithm in the 
present work. In this method, the search direction is 
established as = - fHjJ j where |h/| is the Hessian 

matrix and yf^ Is the gradient vector (~ , ) ^, 

9 2^2 ax 

Along the search direction S,- the optimal step length oc at 

X s 

which the function is minimum in that linear direction is 
calculated using any one of the available one dimensional 
minimization methods, Fibonaeci method is used in this work 


which makes use of the sequence of Fibonacci numbers for placing 
the experiments in the interval of uncertainities. The number 
of experiments is decided as fifteen to get the accuracy of 
0,00103. Having obtained the search direction, and step 
length, oCg , wo move to the new design point according to 

the expression, ~ ^ ^ » where is the old design 

point. The Hessian matrix is updated in each iteration of 


unconstrained minimization as » where , 

[“ij= 


(“.SpCCCgS.) 


{o sT) 

^ S 1 




and 

The flow chart for the SUMT algorithm is shown in Fig. 4. 
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Inorder to have an efficient algorithm for linear 
minimization to obtain, oc* , the initial interval of 
uncertaintly for bracketing the minimum along that direction 
is to be very carefully fixed. The following procedure has 
been used to automatically fix up the initial interval of 
of unC'-rtainity in the present work. 


1, find the values of <x . where, oc, = 

and Qk is the number of negative elements in the vector of 
search direction. These values represent respectively the 
limiting values of step length to avoid the non-negativity of the 
corresponding variable in the next design point. 

2, Select the minimum value among >“ 2 * * • ‘^Ik check 
whether the design point with this minimiun value of this step 
length IS feasible or not. Set this minimum cc as 

3, If the design point is feasible then the initial interval 

of uncprtainity is where 1^=0, 0 and If it is 

not feasible, reduce it to cc = c,a where c is the reduction 
factor and check for the violation of constraints and repeat 
the process till the value a^^jLn obtained which satisfies non- 
negativity of design variables and the constraints. 


4, In some cases there may be a search direction which has 
no negative element. In such a situation any interval of 
uncertainity apriori can be fixed. 
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5» III "the reduction, process, once tlie feasi’ole region 

approaches zero or a very small quanity then a is set to zero. 

The above strategy to select the initial interval of 
uncerbainity is described in the flow chart given in Pig,5» 

4.5 Program Development: 

Iny structural optimization program , particularly 
nonlinear programming problems, basic ally consists of the 
following main parts. 

1, Routine which evaluates the value of the objective function 

2. Routine for carrying out the analysis of the structure 

5. Routine for checking of the constraints, 

4. Routine for the solution technique to reach the optimum. 

In this work, three separate computer programs are 

developed which have similar format except for the difference 
in the corresponding methods of analysis and their objective 
function and constraints. 

following are some situations in which the programs 
need certain modifications, Ihen the balcony slab is provided 
at the ring beam 2 level, the corresponding moment and 
forces are to be considered in the analysis and the necessary 
modifications can he made accordingly. When the cross-sections 
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FIG-5. FLOW CHART TO FIND INITIAL INTERVAL 
OF UNCERTAINITY 
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columns are required as square or rectangular the design vector 
and the necessary changes in the analysis can he easily done. 

For Certain combinations of the number of columns and bracings 
the present analysis of the staging may turn out to give crude 
values. For analysis of the staging with less number of 
columns, say four, and/or for less number of bracings, the 
present analysis should not be used. 

Program 1 : Optimal design of the container: 

This program consists of number of routines which are 
explained below; 

Subroutine MilN: This contains all the necessary data for 

the design and analysis of the container. The values of (1) the 
initial design point, (2) the maxmum and minimum values of the 
design variables for normalization, (3) the cost details of 
concrete, steel and shuttering, (4) the allowable stress values 
of concrete and steel (5) the conical shell coefficients which 
are required to evaluate stiffness (6) other preassigned 
parameters and capacity, number of columns are all given in 
DATA Cards. They are transfered to the required subroutines 
through OOBMON statements, ipart from setting up and 
transferring the necessary values to other routines, this main 
calls SUMT (algorithm) routine with the initial design point 
which returns the optimal design point. 
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Subroutine SIM: This routine forms the unconstrained function 

(0) using the ob;)ecti've function, constraints and the appropriate 
penalty parameter and calls DTP to get a new design point, 
me penalty parameter r is reduced in each iteration and 
the corresponding unconstrained function is minimized using 
PPP. It returws the optimal design point to the main program. 

Subroutine DPP: This routine calls the other routines,DIPP, 

SEAROH, PUNG, PIBBN, HEiSIM. This routine essentially does 
the calculations as shown in Pig*4« "With the starting Hesian 
matrix which is a unit matrix and the gradient vector at 
the starting point (using DIPP)±fc calls SEARCH and obtains the 
Optimum step length from FIBBN, It returns the optimum design 
point to SUMT once the convergence criteria is satisfied. 

Subroutine DIPF: This routine calculates the gradient vector 

using the central difference formula. This calls FURC to 
evaluate the function values at the f 02 ?ward and backward point 
for each design variable. The probe length EPQ is taken as 
1 percent of the design variable value (EPQpO.01) . EPQ 
Can be altered if the higher or lower probe length is derived. 

Subroutine SEARCH: This routine returns the normalized 

search direction vector The Hesian matrix and the gradient 

vector (— yf)are multiplied in this subroutine and is added with 

the previous search direction to get the new search direction. 
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Subroutine HESIAN: This routine does the necessary calculations 

for updating the Hesian matrix. The required quantities for 
updating the matrix are the current Hesian matrix fn.! , the 
difference of gradient vector (Q^), search direction vector (^) 
and the optimum step length (a^). The additional matrices 

generated and added with^Hj^l updated Hessian Matrix 
[HjL_^^]is returned to the routine DFP. 

Subroutine PIBM: This routine returns the optimum step 
* 

length the search direction. This routine calls PHIE 

O 

which evaluates the function value at any point along the search 

direction. This routine implements the strategy to select 

the proper initial interval of uncertainity as explained earlier, 

■X" 

and then searches for oc using the fibonacci search technique. 

Subroutine PHEP: This routine returns the function value ( 0 ) 

at any point ALFA along the search direction S from the given 
point This calls PUNC to evaluate the value and sends 

the value to PIBBN, 

Subroutine PUHC: This routine evaluates the unconstrained 

function value ( 0 ). It calls the routine PN to evaluate 
the objective function value and calls OOFSUB to get the 
summation of tho inverse values of the constraints. The 
necessary penalty parameter is transferred via C0MI40H from 


SUIT. 
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SuTDroutine FH : Ihis routine returns the oh^ective function 
value for any given design point, it calculates the concrete 
and steel quantities and shuttering area of the various elements 
and using the correspoming cost values it evaluates the 
objection function* The necessary invarient parameters are 
transfered through COMiOh statement. 

Subroutine CONSUE: This routine containts the geometric 

constraints and the behaviour constraints and it returns 
the constraint values and the summation of the inverse values 
of the constraints to other routines. If calls ANililS to get 
the values for some quantities in the behaviour constraint 
expressions, GJ(n) represents the nth constraints and GJX 
represents the summation. 

Subroutine iHillS: This routine is developed as a separate 

program for the analysis of the container of the intze tank. 
Both the membrane analysis and continuity analysis are carried 
out in this routine. For a given set of values of the 
variables and other assigned parameters this routine gives 
the forces and moments and the resulting stresses and the 
corresponding reinforcement areas required. The values of 
stresses, steel areas, and volume are returned to CONSUE 


from this routine , 
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Subroutine SOLVE: This routine returns the solution of the 

system of eq.uations giving thereby values of displacement and 
rotation (LISP and ROT) corresponding to the compatibility 
req^uirements. This is called thrice by MXLYS for each set 
of rotation and displacement g,t the three ring beam junction, 

LU factorization technique with pivoting stragegy is implemented, 
in order to avoid numerical instability. 

Subroutine GURBCO: This routine returns the coefficients ^^>^2 
andp^ to calculate the design moments in a circular beam which 
is supported by RSUP or SUPR number of columns. Since the 
coefficients are dependent only on the number of columns which 
is a preassigned parameter , this is called only once by the 
main program and transfered to MAIIS via COMMON, 

Program 2: This program for the optmal design of the staging 

has the samilar routines performing similar operations. The 
main program, analysis routine, routine for calculating the 
value of the objective function and constraints are all 
corresponding to the problem of the supporting frame, whereas 
the routines for SUMT algorithm remains the same. 

Subroutine MAIN: This routine has all the values required for 

the other routines. The initial design point vector, maximum 
and minimum limits for normalization, other preassigned parameter 
like the height of staging, number of columns and panels are 
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all available m. DATA cards and are transfered to the required 
routines via COB/UION statements* The values required for the 
resultant wind force and the distance at which it acts, are 
all given as input which are derived from the optimal design 
of the container (from program 1), 

The other subroutines from SUIT to PHTP are exactly 
the same in this program. 

Subroutine FN: This routine calculates the concrete and steel 
quantities in the columns and bracing beams and then it evaluates 
the value of the objective function using the costs of concrete 
and steel and returns this to FUPC, 

Subroqfcine OOPSUE; This routine evaluates the values of all the 
constraints (given in the earlier chapter) and also the summation 
of the inverse of the constraints. The number of constraints 
depends on the number of panels HP, 

Subroutine ANILYS: This routine performs the analysis of the 
supporting frame both for wind force and earthquake loading. 

The bonding moment and twisting moment which gives maximum 
values for earthquake or wind force are considered as design 
forces’* Similarly the s^Im of stress ratios is calculated for 
both wind and earthquake cases, and the maximum one is returned 


to COWSUE, 
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iTOgram 3: Tins p'^ogram which also has similar routines as 
in the above two programs is to find the optimpm design point 
j'or the foundation. 

MAIN : This has the values of the loads that are to he 
transfered to the soil, the hearing capacity of the soil, and 
the initial design point and the limiting values of the variables 
for normalisation. It sands the values of the curved beam 
coefficients to ANAIYS which are obtained by calling CURBCO, 

Subroutine FN: This calculates the objective function value 

i.e. the cost oi bhe foundation at any design point by evaluating 

the concrete and steel quantities and using the respective cost 
values. 

Subroutine Aj'IAIYS; It performs the analysis of the ring beam 
and raft slab and calculates the required area of steel and 
the corresponding stresses in concrete and steel at the -points of 
maximum force or moment. It requires the curved beam coefficients 
which are transfered from MAIN program. 



CHAPTER 5 


RESULTS, DISCUSSION AND CONCLUSION 
5.1 Introduction: 

Two problems as detailed below respectively of 600 KL 
and 1000 KL cnpacity have been solved to obtaih the minimum cost 
design of R C Intze tank by the use of the developed programs. 
All the results given have been obtained on the IM 7044/I401 
computer at IIT Kanpur, The bearing capacity of the soil 
has been taken to be 20 t/m in each case. The wind pressure 
on the structures is assumed to be 150 Kg/m . ill the allowable 
stresses considered in the present work are as per IS Code, 


5,2 Details of the Problems: 

Problem 1 : Capacity of the tank : 600 KL 
Height of Staging : 15 M 

Concrete Used : M200 

Problem 2: Capacity of the tank : 1000 KL 

Height of Staging : 22 M 

Concrete used for : M200 


Container and Foundation 


Concrete used for 
Staging 


: M2 50 
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Ulie pairoins ui“ic si^udy has hesn. earned oulj hy changing 
"bhe nujiiher of colmnns and "the nuinbei? of panels in "the staging* 

For 600 KL capacity towers, the numher cf columns considered 
in the present work are 6, 8 and 10 and the number of panels 
for its 15 m staging have been taken as 3, 4 and 5 in each 
case. For 1000 K1 capacity tower, the number of columns 
considered are 8, 10, and 12 and the number of panels for its 
22m staging are 4, 5 and 6 in each case. 

For brevity, these problems are refered or designated 
as follows; 

0(600,6) means the container problem of 600 KL capacity 
with 6 columns. 

8(600,6,4) means the staging with 6 columns and 4 Panels 
supporting the 600 KL capacity container. 

F(600,6,4) means the foundation for the 600 KL capacity 
container tank with the staging contains 6 columns and 4 panels. 

5.3 Results: 

The minimum cost of the water tank is obtained using 
the programs 1, 2 and 3 in an interactive manner. That is, for 
a given capacity, the optimal design of the container is 
obtained first, using program 1. Then program 2 is energised where 
the loads coming on each column are calculated using the 
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EABLl 2 t C( 1000, to): OOFmiHER OPTIMIZATICIT KSSUITS 


X 

Starting X" 

^■7 

4 

4 

T 

' 4 ■" 


13.75m 

13.66 

t3.66 

13.65 

13.65 

13.65 

^2 

6 *00111 

5.99 

3.99 

6,00 

6.00 

6.00 

Xj 

35.00 cm 

41.80 

41:.80 

31.50 

31.50 

31.50 

^4 

60.00 cm 

62.50 

62.50 

60*34 

60,33 

60.33 

X5 

85.00 cm 

85.60 

85.^0 

84.16 

84.15 

84.15 


65.00 cm 

69.20 

69.20 

62.87 

62.85 

62.85 

Xrj 

1.70 m 

t.70 

1.703 

1 ,7P8 

1.707 

1 .707 

Xq 

10,00 m 

10.00 

10.00 

10.03 

10,03 

10.03 

2C9 

60*00 cm 

49.20 

4.9.2 

46.5 

46*52 

46,52 

^10 

100,00 cm 

98.80 

95.80 

97.86 

97.96 

97.96 

^1.1 

20,00 cm 

20.90 

20.90 

14.57 

14.57 

14.57 

^12 

0.004 

0.021 

0*0968 

0.0033 

0.0031 

0.0031 

^13 

0.004 

0.014 

0,004 

0.0032 

0.003 

0.003 


0.004 

0.018 

0,003.8 

0.0032 

0.0031 

0.0031 

^15 

0,01 

0.004 

0.0039 

0.0037 

0,0031 

0.0031 

V 

“ 

402 

40.2 

4.02 

0.4 

0.04 

%; 

118316*67 

165227.8 

1222802.3 

104592.1 

103705*1 

103705.1 

1,11., 12, 15, ■ 

'8, It, 14,1 

8,12,13, 

... 8,-1t,12 

, 12,8, 

'12,8, 11, 

Ac tive 

22,23,24 

3t 

14., 1, 

13,14, 

11,13,14 

13, 14,22 

Constraints 3t 


31 

22,23, 

22,23,24 

23,24,1,: 




ziai,. 
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TABLE 3 s PEOBLEl S(600,TO,4) s OPTIldZATICM OHl STAGIBO 


Design 
Vector X 

Starting iT" 


4' ' 


4 


50*00 cm 

49*89 

4,9.89 

49.8 

49.8 

^2 

25.00 om 

24.90 

24.90 

24.90 

24.9 


60.00 cm 

59.90 

59.90 

59.90 

59.90 

^4, 

0.01 

0,015 

0.00818 

0.0089 

0.0086 


0.01 

0,0085 

0.035 

0.0085 

0,0081. 

^6 

0.01; 

0.014 

0.0108 

0.00812 

0.0084, 

Xj 

0.01 

0,015 

0.0818 

0,0818 

0,008 

X 

8 

0.015 

0.0053 

0.0037 

0.00313 

0.003 

Xg 

0.015 

0.0053 

0.0035 

0,00307 

0,00314 

^1.0 

0,015 

0.0053 

0,0034 

0.0033 

0.00312 

r 

- 

64*37 

6.43 

0.643 

0.064 

Pj^(Es.) 25.35 1..33 

23298.53 23821.08 

18965.91 

18875.59 


Active 

Const- 

raints 


5, 6, 8, 9, 9,5, 6, 8, 9, 5, 6, 8, 1.0, 9,5, 10,1 1., 9,5,1.0,11,12., 

10,11, 10,11 11,12,1:3, 1i2,13>14, 13,14,6,8 

14 6,8 


2 >- penalty parameter 
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TABLE 4: PRCBIEM P(600, ) ; POMLATOT OPTIMIZATIOIT 


Besiga 

Vector 

Starting X 





52*0 cm 

55.67 

53*4 

53.4 

^2 

90.0 cm 

54.49 

57.09 

55.09 


25 *0 cm 

29.53 

29.51 

29.50 

^4 

20.0 cm 

tf7.08 

17.03 

17.00 


100.0 cm 

94.92 

94.32 

94.28 

^6 

0.01 

0.0161 

0.011 

0.0109 


0.005 

0,0103 

0.0055 

0.0053 


0.01 

0.0174 

0.015 

0.004 

Xg 

0.008 

0.0065 

0.0063 

0.0062 

r 

.. 

19.48 

1.948 

0.194 

) 

12455.58 

10941.08 9444.08 

9420.88 


1y2,7,15, t, 2,7, 11, 1,2,7,11, 

1,8,19 14,15,16, 14,15,16,18, 

18,19 19 


Active 

Const- 

raints 


1,6,7,18, 

19 
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TABLB 5 i GCMTAISm OPTIlvUZATIOU SESULTS 


Design 

Vector 

— ijfc. 

0(600,6) 

0(600,8) 

0(600,10) 

Starting 

Point 

Optimum 

Point 

Starting 

Pomt 

Optimum 

Point 

Starting 

Point 

Optiitum 

Point 

X^(m) 

12.25 

12.18 

12.25 

12.24 

1i2.25 

12.24 


4*50 

4.47 

4.50 

4.43 

4.50 

4.43 

(cm) 

25.00 

27.52 

25.00 

28.09 

25.00 

28.07 

X^(cm) 

40.00 

29.59 

40.00 

38.89 

40.00 

38.80 

(cm) 

45.00 

56.80 

45.00 

51.63 

45.00 

51.60 

Xg (cm) 

45.00 

45.15 

45.00 

47.96 

45.00 

47.90 

Jrj (m) 

1.50 

1.42 

1.50 

1.426 

1.50 

1.426 

X3(m) 

9.00 

8,96 

9.00 

8.975 

9.00 

8,98 

Xg(cm) 

45.00 

20,02 

45.00 

20.21, 

45.00 

20.20 

^ho(cm) 

80.00 

76.99 

80.00 

74.28 

80.00 

74.28 

X^ ^ (cm) 

20.00 

15.06 

20.00 

11.83' 

20.00 

11.80 

^12 

0.004 

0.0152 

0.004 

0.0032 

0.004 

0.0036 

X^3 

0.004 

0.0170 

0,004 

0,0058 

0,004 

0.0038 

^14 

0.004 

0.004 

0.004 

0.003 

0.004 

0.003 

^15 

0.02 

0.0072 

0.01 

0.004 

0,01 

0.0038 

^1 

350.00 

- 

273.00 

~ 

273.80 

- 


~ 

0.350 

- 

0.273 

- 

0.273 

%(Rs.)7 1356.7 

70721.09 

68227.78 62298.67 68227*78 

62329.67 





OJable 5 

c ontd . * 4 
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Table 5 contd»**» 


Design 
Tec tor 

— ‘V 

X 

- 

0(1000, 6)) 

0(1000,12) 

Starting 

Point 

Optimum 

Point 

Starting 

Point 

Optimum 

Point 

Z^(m) 

13.75 

13.65 

13.75 

13.68 

ZgCni) 

6.00 

6.00 

6.00 

5.97 

X^Com) 

35.00 

31.52 

35.00 

43.28 

I^(cni) 

60,00 

60.33 

60.00 

62.12 

I^(cm) 

85.00 

84.15 

85.00 

85.50 

Xg(om) 

65.00 

62.85 

65.00 

69.74 

S^Cm) 

1.70 

1.707 

1.70 

1.69 

XgCm) 

10.00 

10.03 

10.00 

9.98 

Xg(cm) 

60,00 

46.52 

60.00 

3.5.66 

X^Q(om) 

100.00 

97.86 

100.00 

96.50 

X^^(cm) 

zo.oo 

14.54 

20.00 

23.39 

^12 

0.004 

0.01,6 

0.004 

0.0039 

Xi3 

0.004 

0,0169 

0,004 

0.00301 

^14 

0.004 

0,0039 

0.004 

0.00306 

^15 

0.02 

0,0076 

0.01 

0.0031 


402.00 

- 

282.00 

- 

* 

r 

~ 

0.40 

- 

0.28 

Puncti on 
Values 
(fe. ) 

126058.52 

123987.60 

118265.32 

117822.00 


- Initial penalty parameter; r^-pmal penalty parameter 
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TiBLB 6 : STAGING OPriMZATION RESULTS 


Design 

Vector 

3(600,6,3) 

S(600,6,4) 

S(600, 8,4) 


Starting 

Point 

Optimim 

Point 

Starting 

Point 

Optimnm Starting 
Point Point 

Optimum 

Point 

X^(om) 

60.00 

63.45 

60.00 

64.36 

55.00 

53.30 

XgCom) 

0 

0 

C\j 

24.45 

25.00 

24.60 

25.00 

24.39 

X^(cm) 

60,00 

41>.90 

60.00 

59.80 

60.00 

59.70 

24. 

0.02 

0.00838 

0.02 

0.0082 

0.02 

0.0079 


0.02 

0.0087 

0.02 

0,00827 

0.02 

0,0087 

% 

0.02 

0.00802 

0.02 

0.00849 

0.02 

0.0081 

Zj 

0.015 

0,0039 

0.02 

0.008 

0.02 

0.009 


0.015 

0.0030 

0.15 

0.00303 

0.015 

0,0031 


- 

- 

0.015 

0.00303 

0.015 

0.0031 

ho 

- 

- 

0.0t5 

0.00302 

0.015 

0.0032 

■^^11 

- 

- 

- 

~ 

~ 

- 

^12 

— 

— 

- 

- 

— 

- 


344.1 

- 

210.00 


347»60 

- 

* 

r 

- 

0.344 

- 

0.210 

- 

0.347 

Pij(Rs.: 

)25624.3 1 6054 . 8 

28903.0 

18535.90 31169.6 

17630.4 


QlaUle 6 contd 
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Table 6 contd 




Design 

Vector 

> 

' S(600,8,5) ■' 

3(600,10,3) ' 

' 3(600,10,5) 

S^-taarting 

Point 

Optimum 

Point 

Starting 

Point 

Optimum 

Point 

Starting 

Point 

'Optimum 

Point 

(cm) 

55.00 

54.98 

55.00 

47.90 

55.00 

54.60 

XgCcm) 

25.00 

Z4.89 

25.00 

20.32 

25-00 

24.80 

(cm) 

60.00 

59.90 

60.00 

57.96 

57-96 

59.90 

^4 

0.02 

0.008 

0.02 

0.00804 

0.02 

0.008 

^5 

0.02 

0.0085 

0.02 

0.0081 

0,02 

0.008' 


0.02 

0.008 

0.02 

0.0087 

0.02 

0,008( 


0.02 

0,0082 

0.015 

0.0034 

0.02 

0,008 


0.02 

0.0080 

0.015 

0.0032 

0.02 

0.008'! 

Zg 

0.015 

0.0031 

- 

- 

0.015 

0.003 

^10 

0.015 

0.0035 

- 

- 

0.015 

0.032 

^11 

0.015 

0.0035 

~ 

- 

0.015 

0.005c 

^12 

0.015 

0.0036 

- 


0,015 

0.0035 


157.9 

- 

630-9 

- 

437.50 

- 

¥: 

X 

- 

0.157 

- 

0.0636 

- 

0.0436 

Pu(Bs.)34414-0 

19930 .0 

33046.52 

15820.51 39298.20 23247.69 
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TABIjE 1 : STAGIIG OPTIMIZATION: EESUITS 


Design 

Vector 

X 

S( TOGO, 8, 4.) 

S(T000, 

8,6) 

3(1000, 

10,6) 

Starting 

Point 

Optimum 

Point 

'Starting 

Point 

'Optimiim 

Point 

St^'tirig' 

Point 

Optimum 

Point 

(om) 

60.00 

54.80 

60,00 

53.50 

60.00 

56.00 

XgCom) 

15.00 

21.20 

25-00 

21 w 50 

25.00 

23.70 

X^ (cm) 

60.00 

57.60 

60.00 

58.38 

60.00 

59.40 

^4 

0.03 

0.03 

0.03, 

0,0083; 

0.03, 

0.0084 

Xe 

5 

0.03 

0.008 

0.03; 

0.057 

0.03 

0.016 

^6 

0.03 

0.008 

0.03, 

0.0086 

0.03; 

0.011 


0.03 

0.0084 

0.03; 

0.0082 

0.03 

0.0088 


0.015 

0.0033 

0,03, 

0.008 

0,03 

0.0083 

Xg 

0.015 

0.003 

0.03 

0.0081 

O.O3 

0.0084 

^10 

0.015 

0.003, 

0.015 

0.0031 

0.015 

0.0035 

^11 

- 

- 

0.015 

0.013 

0.015 

0.003' 

Xi2 

- 

- 

0.015 

0.015 

0.15 

0.0031 


- 


0.015 

0,0035 

0.015 

0.003 

Zl 4 

— 

- 

0.015 

0.010 

0.015 

0.0033 

^1 

108.2 

- 

79.2 

~ 

102.0 


* 

r 


0.01 

- 

0.079 

- 

0.01 


) 60415.6 

30820.0 

67878.92 39273.81 

79528.1 

38443.97 


Table 7 contd 
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Design 
Vec ’cor 

3(1000, 12,4) 

S( 1000, 12,5) 

Starting 

Point 

Optimim 

Point 

Starting 

Point 

Optimiim 

Point 

(om) 

60.00 

58.00 

6O0OO 

59.00 

X^Ccm) 

25.00 

24.70 

25.00 

24.80 

Z^ (cm) 

60.00 

59.50 

60.00 

59.90 


0.03 

0.0081, 

0.03 

0.0085 


0.03 

0.009 

0.03 

0,01 

"^6 

0,03 

0.009 

0.05 

0.008 

X, 

0.05 

0.008 

0,03 

0.0083 

Zg 

0.015 

0,0048 

0.03 

0.00817 

Xg 

0,015 

0,0044 

0,015 

0.0138 

2=10 

0,015 

0,02 

0,015 

0,012 

tr 

"~11 

- 

- 

0,015 

0,0056 

^12 

~ 

- 

0.015 

0,0065 

Xi 3 


- 


- 

2-14 

— 

— 

— 


^1 

165.70 

- 

159.0 

- 

r* 

- 

0.165 

- 

0.139 

84182.17 

44976.1 

87632.5 

46320.40) 
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FOUNDATION OPTir/CEZAHON EESIILTS 


Design 

Vector 

— 

X 

fC^oo, 6,3) 

P(600,6,4,) 

P(600,8,5) 

starting 

Point 

■ Optimism 
Point 

Starting 

Point 

Optimum 

Point 

Starting 

Point 

Optimum 

Point 

(om) 

65.00 

77.40 

65.00 

69.32 

55.00 

55.90 

XgCcni) 

100.00 

93.56 

100.00 

77.88 

90,00 

70.40 

(c]^) 

30,00 

26,96 

30.00 

38.18 

25.00 

20,00 

Z^(cm) 

25.00 

16.01 

25.00 

15.00 

20.00 

14.99 

Z^Com) 

150,00 

143.00 

150.00 

14f5.05 

100.00 

84. 6^ 

Tr 

^6 

0.02 

0,0097 0.02 

0.011 

0.01 

0.0032 


0.01 

0.0049 0.01 

0.0056 

0.005 

0,003 

% 

0.02 

0.013 

0.01 

0.0043 

0.01 

0.0166 

Zg 

0.01 

0.0084 0.01 

0.0062 

0.008 

0.003 

^1 

44.80 

- 

42.4 

- 

21.59 

- 


— 

0.44 

- 

0.004 

- 

0.002 

1 1 

C- 

• 

26259.2 

18297.17 

26255.1 

16987.0 

12808.35 

6232.1 


Table 8 contd**..* 
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Design 

Vector 

— -ai*’ 

X 

P( 1000, 8/6) 

■ 

o 

o 

10,6) 

P(1000, 12, ■'6) 

Starting 

Point 

Optimm 

Point 

starting 

Point 

Optimum 

Point 

Starting 

Point 

Optimum 

Point 

(cm) 

55.0 

63.90 

58.00 

67.90 

62.00 , 

67.36 

ZgCom) 

100.0 

92.08 

90.00 

84.63 

90.00 

82.98 

(om) 

50.0 

38.90 

50.00 

43.36 

50.00 

36.38 

Z. (om) 

4- 

:»5.o 

18.51 

25.00 

17. Z8 

25.00 

19.66 

Z^ (om) 

2.00.0 

152,48 

ZOO. 00 

180.94, 

200.00 

156.93 

^6 

0.02 

0.01, 

0.02 

0.0086 

0.02 

0.0075 


0,01 

0.0052 

0.01 

0.00434 

0.01 

0.0037 

Zg 

0.022 

0.0093^ 

0.022- 

0.005 

0.022 

0.01 

^9 

0.01 

0.008^ 

0.01 

0.008 

0.01 

0.009 


65. 8Z 

- 

64.07 

- 

66.37 

- 

r* 

- 

0.658 

- 

0.064 


0,066 

%(^’^f ) 

37417.69 

22070.0 

36579.5 Z4092.6 

37424.0 

20624.4 



80 


■fehe geometry of the optimised container, and wind pressure and 
earthquake requironents are given as data. This program then 
obtains the optimal design of the staging • Pinally^the data fro 
the program 2, like the load coming on the foundation, the 
diameter of the bottom ring beam etc, are used in program 5 
so as to obtain the optimal foundation design, 

5,4 Discussion: 

The pT^oblems which have been solved in che proceeding 
section show the flexibility of the developed programs which can 
be used for any capacity, any height of staging with a given 
number of columns and panels to obtain the optimal design. 

The starting point of the problem 0(600,8) (i.e,) contains] 
of 600 KD capacity supported on 8 columns and the optimal design 
point of the same problem are given in the fable 5 • The 

cost of the container is minimized from Rs, 68227 to Rs. 62298 
which which shows 8,69 percent reduction. The starting point 
for this Case is the design vector worked out in Ref. (8). 

The minimization history shown for the problem 
0(1000,10) shows that the geometric variables (i.e,) diameter 
of the cylinderical wall, , height of the cylinderical 
wall, X 2 y height of the conical shell, X^ , and the diameter of 
the ring bcajii 3 > nre not significantly different from 
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the Starting point and the optimum point. This is because of 
the constraint GJ(8). Since this constraint is a function 
of the variables which affects the volume of the container, 
any small change in these variables which may reduce the 
function value is prevented in order to satisfy this constraint. 
This observation leads us to conclude that the size of the 
problem can bo further reduced by fixing the configuration 
(neglecting as design variables) of the container 

in which case the shuttering cost shall be constant and need 
not be considered as a part of the optimal design formulation. 

It can be seen that the optimal cost of the container 
for a given capacity varies with the variation of the number 
of columns in bhe staging. The optimal cost of the container 
for 600 EX capacity is more for 6 number of columns as compared 
to that of 8 and 10 columns. Similar observation is made for 
1000 KL capacity also. That is, the cost for 8 columns is 
more as compared to that of 10 and 12 columns. This is because 
the bending moment in the ring beam 3 is more when the number 
of columns is less. 

The optimal design of the staging for problems are 
shovm in the Table 7 , It can be seen that the minimization 

is achieved by the reduction of the steel in the columns and 
bracing beams towards the minimum reinforcement requirement 
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as per IS code and diameter of the column is getting fixed 
accordingljr. In the staging where the number of columns 
are lesser say 3(600,6,4) , the diameter of the column 
becomes larger at the optimum design corresponding to minimum 
reinforcement, in order to carry the load coming on it. It 
IS obvious that the load coming on a column is reduced when 
the staging contains more number of columns. However, it can 
be seen from the Table 7 that the minimum cost of the staging 
are not significantly different as the number of columns vary. 

But the cost of the foundation varies significantly for the 
variation in the number of columns because the cost of the 
foundation depends on the load coming on the column(load per 
span on the bottom ring beam between the two columns) and also 
on bhe self weight of the staging. For 600 KL capacity problem 
the cost of staging for 3(600,6,4) , 3(600,8,5) and 3(600,10,4) R* 
are around Rs. 19,000 but the foundation cost of the corres- 
ponding problons show significant variation for the obvious 
reason explained above. Thus it can be concluded that the 
optimum design of the staging can be obtained when the 
staging is designed with more number of columns i,e. 10 for 
600 KL capacity and 12 for 1000 KL capacity. 

The solution of the staging for problem 2 are shown 
in the Tablo 7 . 3ince 60 cm column with 3 percent of steel 
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was found unsafe with M200 mix , design with M250 is attempted* 
It violates the constraint GJ(4) for the stress limits, that is, 




< 


( Refer . 


Inst-ad of trying with larger diameter with the M200 mix, 
the design using M250 concrete mix is attempted. In the case 
of the staging for larger capacity containers it is better to 
go in for the richer mix for staging because it may reduce the 
weight of the staging which inturn will reduce the cost of 
the foundation. It is obvious from the solutions of 600 KL 
capacity problem, that the diameter of the column will become 
higher with the steel towards the minimum which increases the 
weight of the staging though its cost is reduced. It can be 
concluded that the staging design should be attempted for 
different mixes in order to pick up the better minimum design 
of the entire structure, keeping in mind that the foundation 
cost depends on the weight of the staging. 


The solution of the foundation design problems are 
shown in the Table 8 . It can be seen that the minimization 
of the cost is achieved by the reduction of the slab width, 3^, 
reduction of the thickness of slab, towards the minimum 
thickness reduction in the depth »f the ring beam, Xg, and the 
reduction in steel variables. In problems I'(600,6,4) and 
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P(600,'8»,4) the thickness of the slat at the junction of the 
heam, X^, is increaded from the starting point value. This is 
because of the increase in soil pressure due to the reduction 
in the hearing area as the slab width decreases. 

The variables and are the steel variables expressed 
as fraction of the cross section. It can be seen that in the 
solutions at the optimum point, the values are in the ratio of 
2:1 approximately, though they are treated as two separate 
variables in the design vector. This behaviour is one which 
is expected because the negative and postive bending moments 
are also in the same ratio. In the formulation of the container 
problem one of these two steel variables (in ring beam 3) only 
is included in the design vector unlike both in the foundation 
problem. Thus the behaviour of the variables in the solution 
of the foundation ring beam assures the correctness of the 
formulation of the container with one steel variable in the 
ring beam 5. 

It has been observed from the optimised values of the conts 
iner that the total loa^ transfeced to the staging^ and the 
diot.mco at which equivalent earthquake force o^ win^orce 
on the coubaineji^are not 4Ur ffe ' 3 ? in § significantly^^,!^ different 
number of columns or braces, Prom this it can be 

concluded that there is no interaction between the container 
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and ^ staging, “ bU ' t iti ' in there -^hetween staging and foundation. 

But in the present work the staging and the foundation are 

seperately optimised and to that ejcrtent the results obtained 

are approximate. 

Cost Vectors and Cost Matrices 
1, Container cost vector of 600 KL Problem: 


Columns Cost in Rs, 

6 70721 

8 62298 

10 62329 


2. Container cost vector of 1000 KL problem: 


Columns Cost in Rs, 

123,987 
1,03,705 
1 , 17,822 


8 

10 

12 


The cost values in the following matrices are in rupees 
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3* Staging cost matrix of 600 KE Problem: 


Panels 

3 

4 

5 

Columns 

6 

16,054 

18,535 

20,150 

8 

15,990 

17,630 

19,930 

10 

15,820 

18,875 

23,247 


4-. Foundation cost matrix for 600 EX problem: 


Panolg 

Oolumns 

3 

4 

5 

6 

18,297 

16,987 

17,090 

8 

10,459 

9,420 

6,232 

10 

12,320 

11,340 

10,200 


5« Staging cost matrix for 1000 KX problem 


Panels 

Oolumns 

4 

5 

6 

8 

30,820 

34,450 

39,276 

10 

35,940 

37,120 

38,443 


44,9 76 46,320 49,300 


12 



6« Foimdation cost matrix for 1000 KL problem* 


Panels 

4 

5 

6 

Golumns 


8 

38,229 

37,625 

22,070 

10 

37,245 

36,725 

24,092 

12 

21,941 

20,895 

20, 624 


Total cost liatrix for 

600 KL problems 


Panels 


' 

— 3 

4 

5 

Columns 




6 1 
8 

TO 


,05,072. 1,06, 
88,747 89, 

90,469 92', 


243 

1,09,961 

348 

88,460 


(minimum) 

544 

95,776 


Total cost matrix for 1000 KL problems 


Panels 4 5-6 

Columns 

8 1,93,036 1,96,062 1,85', 333 

10 1,76,890 1,77,550 1,66,240 

(minimm) 

1,84„729 1,85,037 1,87,746 


12 
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Prom the total cost matrices it can be observed that the 

minimum cost design for 600 Kl problem is the one in which the 

! 

staging contains 8 columns and 5 panels. Similarly, the 
minimum cost design for 1000 H, capacity problem is the 
one in which the staging contains 10 columns and six panels. 

That IS, the optimum design of the entire structure for 600 KL 
capacity problem can be had from 0(600,8), 3(600,8,5) and 
P(600,8,5) , and for 1000 KL capacity problem the optimum 
design can be had from 0(1000,10), 3(1000,10,6) and P( 1000 , 10,6) . 

5.5 Conclusions: 

Prom the above discussions, the following conclusions 
are summarised, 

1 , The geometric variables in the container are not 
varying significantly to obtain the minimum cost. 

2, The reinforcements in the container ring beams and 
the walls , and the reinforcement in the columns go to the 
giinimum at the optimum design point, 

3, Richer mix in the staging for large capacity container 
gives the better design. 

4.. Por 600 KL capacity problem the design with 8 columns 

and 5 panels and for 1000 KL capacity problem the design with 
10 columns ond 6 panels are found to give the minimum cost designs 
5, The interaction between the container and the staging 

is obcoi'vcd to be less significant. 
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5.6 Proposed Extension of the fork: 

0}he solution of any investigation ends with certain 
Imitations and apprccximations. But the solution ?;ith 
acceptable mismatch from the exact solution can measure the 
efficiency of the investigation. In this work, the approximatio 
and the limitations of the methods and the formulations are 
pointed out then and there. Pollov/ing are the problems which 
Can be investigated to get better study over this work, in 
a better computing environment. 

1. The separate programs which are used to get the optimal 
design of entire structure can be coupled together to get the 
solutions in automate manner with the provision of selecting 
the proper type of foundation and use the same in the 
corresponding routine, 

2. The space frame program can be developed and the 
approximate analysis of the supporting frame with the exact 
solution Can be comparad. The validity of the approximate 
analysis of the frame, with the variation of the number of 
columns and number of bracing s can be investigated. 

3. In the optimum design of the staging the panel heights 
are kept equal. The optimum cost of the staging with unequal 
panel heights c,m be investigated keeping the panel heights 
as dcsif-gn variables. 



4. The numher of columns and the number of bracings can be 
considered as integer variables and the solution can be 
sought using mixed integer programming algorithms. 

5. Rather than carrying out substructure optimization 
and solving each one of them as RIP, the entire problem can 
be handled as a 3 stage Dynamic programming problem. 
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Pollowajig are th.e cost values used in tins present worlc: 


Oost of concrete in the container (J4200) 
Cost of concrete in the staging (ivl^OO) 
Cost of concrete in the staging (15250) 
Cost of concrete in the foundation (11200) 
Cost of steel 

Cost of straight shuttering 
Cost of conical shell shuttering 
Cost of spherical shell shuttering 


: Rs, 3OO per 

: Es* 280 per nP 

t Es. JOO per nP 

i Es. 260 per mr 

i Es. 2800 per tonne 
t Rs. 20 per 

2 

s Es. 25 per m 
: Es. 50 per 



